On a dual model for 5/ electrons 
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We study the effect of intra-atomic correlations on anisotropics of hopping of 5/ electrons. It is 
shown that they may considerably enhance these anisotropics to the extend that electrons in some of 
the 5/ orbit als remain localized. This provide a microscopic basis for a previously made assumption 
that some of the 5/ electrons must be partially localized when Fermi surfaces and effective masses 
of the component are calculated. Calculations are performed on two- and three-sites clusters. They 
include phase diagrams in the presence and absence of an external magnetic field as function of 
anisotropic hopping. 
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I. INTRODUCTION 



Systems with strongly correlated electrons are a sub- 
ject of considerable research efforts. Among those are the 
ones with 5/ electrons, i. e., involving actinide ions of 
particular interest. While 4/ electrons, because of their 
closeness to the nuclei, are correlated as strongly as pos- 
sible on an atom, the 5/ electrons have a larger spatial 
extent and tend more towards derealization. The de- 
gree of derealization is even more pronounced when 3d 
electrons are considered. By continuing this line of argu- 
ments one finds that s electrons are the least correlated 
ones on an atom and have the largest tendency to delo- 
calize. 

Returning to 5/ electrons, the question arises under 
which circumstances localized features dominate delocal- 
ized ones or vice versa. In particular we want to follow 
up the suggestion which has been made previously [1] 
that 5/ electrons have a dual character in the sense that 
some of the 5/ orbitals lead to derealization when oc- 
cupied by electrons while others result in localization. It 
was demonstrated that a division of 5/ orbitals into de- 
localized and localized ones is able to explain the heavy 
quasiparticles observed in UPt3 and UPd2Al3 including 
their mass anisotropics [1,2]. However, the question re- 
mained unanswered as regards the microscopic origin of 
this division of orbitals. Apparently, strong intra-atomic 
correlations seem to enhance anisotropics in the kinetic 
energy matrix elements which reflect anisotropics in the 
hybridization with the chemical environment. 

The purpose of the present work is to clarify the phys- 
ical origin of the dual character of 5/ electrons. For that 
purpose clusters of two and three sites are considered. 
Because of the large spin-orbit coupling each 5/ electron 
is characterized by its total angular momentum j = 5/2 
and azimuthal quantum number j z . We do not include 
j = 7/2 states in our analysis since their energy is higher 



by approximately leV [3]. We shall treat clusters where 
the average number of 5/ electrons per site is either 2.5 
or 2.67, values which are close to the 5/ occupancy of U 
ions in compounds like UPt3 [3], UPd2Al3 or UGe2 [4,5]. 

Depending on j z different hopping matrix elements are 
assumed for the / electrons in the cluster. The key re- 
sult of the present investigation is to demonstrate that 
indeed anisotropics in the hopping are strongly enhanced 
by intra-atomic correlations. 

The paper is organized as follows. After defining the 
model Hamiltonian a two-site and a three-site cluster are 
investigated. Various properties are calculated, like the 
ground-state energy and magnetization as well as the 
density of states and its magnetic field dependence. This 
is done for different ratios, of hopping. Partial localiza- 
tion is the subject of the subsequent Section V. A discus- 
sion and the conclusions are found in Sect. VI. The Ap- 
pendices are devoted to perturbation theory, the atomic 
limit of vanishing inter-atomic electron transfer and to a 
more detailed analysis of the ground state's composition. 



II. MODEL HAMILTONIAN 

We investigate a microscopic model of 5/ electrons 
which includes local Coulomb interactions and hopping 
terms between nearest neighbor sites. The model Hamil- 
tonian is 

h = - l u (4 ( n ) c 3* ( m ) + h - c ) 

(nm),j z 



^2 hn i zC \ z ( n ) c i* ( n ) + H ° oul 



(1) 



where the first sum is over neighboring sites (nm). Fur- 
thermore Cj z (n) (cj z (n)), creates (annihilates) an electron 
at site n in the 5/ j = 5/2 state with j z = —5/2, ... ,5/2. 
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We will consider two and three sites models. The effec- 
tive hopping between sites results from the hybridization 
of the 5/ states with the orbitals of the ligands and de- 
pends generally on the crystal structure. Rather than 
trying to exhaust all possible different lattice symmetries, 
we shall concentrate here on the special case that hop- 
ping conserves j z . While this is certainly an idealization, 
it allows us to concentrate on our main interest, i. e., 
a study of the influence of atomic correlations on the 
renormalization of hybridization matrix elements. The 
parameters tj z (= t-j z ) are chosen in accordance with 
density-functional calculations for bulk material which 
use jj z basis states. We also include local magnetic fields 
h(n) to account for the influence of ferromagnetic (FM) 
or antiferromagnetic (AFM) internal fields. The local 
Coulomb interactions can be written in the form 

Uj zl j z2 j z3 j z4 c] zi (n) c] z2 (n) c jz3 (n) c jz4 (n) (2) 
with Coulomb matrix elements 

Uj zl j z2 ,j z3 jz4 ~ 3j zl +j z2 ,jz3+jz4 X 

Y.^li^i^\JJ z )Uj{JJ z \l \UU)- (3) 

J 

Here J denotes the total angular momentum of two elec- 
trons and J z = j z i+j Z 2 = jz3+jz4- The sum is restricted 
by the antisymmetry of the Clebsch-Gordan coefficients 
(f f jzi jz2\JJz) to even values J = 0, 2, 4. We use in the 
actual calculations Uj values which are determined from 
LDA wavefunctions for UPt 3 [1], i. e., U J=4 = 17.21eV, 
Uj= 2 = 18.28eV, and U J=0 = 21.00eV. We expect 
Uj=4 < Uj=2 < Uj=o to always hold for Coulomb in- 
teractions, independently of the chemical environment. 
In contrast, the relative order of the hopping matrix ele- 
ments will vary strongly from one compound to the next. 

The average Coulomb repulsion of about 20 eV is irrel- 
evant for the low-energy physics of the model. It simply 
restricts the relevant configurations to states so that the 
sites are occupied either by 2 or 3 electrons. The low- 
energy sector is solely determined by the differences of 
the Uj values, which are of the order of leV and thus 
slightly larger than typical bare / bandwidths. The lat- 
ter are obtained, e.g., from LDA calculations for metallic 
uranium compounds like UPt3. Note that restricting the 
model to f 2 and f 3 configurations is equivalent to let the 
different Uj — > oo while their differences remain finite. 

To mimic the situation in the U-based heavy-fermion 
compounds we consider the intermediate valence regime. 
Note that in the absence of a magnetic field all states of 
the two-site model with 5 electrons will be at least doubly 
degenerate because of Kramers' degeneracy. 

The Hamiltonian Eq. (1) conserves J z = ^2 n J z (n) 
where J z is the z-component of the total angular mo- 



mentum of the system and the J z (n) refer to angular mo- 
mentum projections on individual sites. We shall there- 
fore characterize the eigenstates by their J z value. For 
isotropic hopping tj z = t = const the system is rota- 
tionally invariant. Then J 2 provides an additional good 
quantum number. 



III. TWO SITE MODEL 



As mentioned before it should suffice to limit ourselves 
to the subspace of f 2 and f 3 configurations, because of 
the large values of Uj. To assess the validity of this sim- 
plification we shall compare the predicted ground states 
for the f 2 -f 3 restricted model with those obtained by 
diagonalizing the Hamiltonian Eq. (1) in the full space. 
The rational for the restriction to f 2 -f 3 is not the reduc- 
tion of the numerical diagonalization effort by a factor of 
about two, but the conceptual simplicity gained for the 
interpretation and presentation of the results for, e.g., 
local spin configurations and orbital occupations. In the 
full model, interactions of the order of P /U are present 
which arise from virtual f 2 -f 3 <-> J 1 -/ 4 transitions (t and 
U are representative values of tj z and Uj, respectively.) 
The restricted model ceases to be a good approximation 
when these interactions become comparable to the bare 
parameters of the model. 

For most of the calculations below we shall assume 
^5/2 = ti/2 an d £3/2 as independent parameters, a choice 
motivated by LDA calculations for UPt3 [1] where t 3 / 2 > 
^1/2? ^5/2- However, partial localization, which we will 
discuss in Section V, is also found for other parameter 
choices, e. g., by choosing t 5 / 2 = £3/2 an d t 1 / 2 as inde- 
pendent parameters. The origin is qualitatively identical 
for both choices, as we will demonstrate later. The results 
presented below are obtained from numerical studies of 
the Hamiltonian (1). 



A. Ground-state energy 



Figure la shows the variation of the ground- state en- 
ergy Eq as a function of the hopping parameters t 3 / 2 
and £5/2 = £1/2 • This energy is smooth except for a kink 
along the isotropic line ti/2 = £3/2 = £5/2 • 
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FIG. 1. (a) Total ground-state energy as a function of £3/2 
and £5/2 =ti/2- (b) Expectation value of the Coulomb inter- 
action (^ gs |ifcoui|^gs) in the ground state ^ gs vs. hopping 
parameters £ 3 / 2 and £ 5 / 2 = £1/2- 

Level crossings are apparent from discontinuous jumps 
of the expectation value Ec of the Coulomb interaction 
with respect to the ground state, see Fig. lb. The 
ground- state energy at vanishing hopping is 69.53 eV 
which is the sum of the Coulomb energy of the J = 4 
f 2 configuration and the J = 9/2 f 3 configuration, 



E G ({t Jz = 0}) = E c (f\ J = 4) + E c (f\ J = 9/2) 
5 \ „ 9 



14 



Uj-_ 



14 



J=2 



(4) 



In general, the ground-state energy decreases with in- 
creasing hopping due to a lowering of the kinetic energy. 
This is however partially compensated for by an increase 
in the interaction energy. The jumps indicate quantum 
phase transitions, which result from the competition be- 
tween Coulomb interactions and anisotropic hopping or 
hybridization. 



B. Magnetization 

In order to characterize the ground state of the two- 
site cluster for different choices of the hopping matrix 
elements tj z we study the total magnetization character- 
ized by J z in the presence of a vanishingly small field. 
Results are shown in Fig. 2. The different ground states 
are separated by solid lines when we allow for f 2 and 
f 3 configurations only. Dashed lines show the ground- 
state boundaries when this restriction is lifted. The later 



are obtained from diagonalization of H in the full space 
of 5-electron states. Because of the small differences we 
shall limit ourselves to f 2 -f 3 configurations only. Five 
different phases are visible in Fig. 2, i.e., two "high-spin" 
phases (J z = 15/2, J z = 11/2) and three "low-spin" 
phases (J z = 5/2,3/2,1/2). In the "high-spin" phases 
which prevail for £ 3 / 2 ^> ^1/2 = ^5/2 the two sites are 
ferromagnetically correlated, (3(a) • 3(b)) > 0, while in 
the low-spin phases the correlations are antiferromag- 
netic (AFM), (3(a) -3(b)) < 0). Here and below we 
denote the two different sites by a and b. 




0.6 0.8 1.0 
t 1/2 =t 5/2 ( eV ) 

FIG. 2. Phase diagram in the £i/2(=^5/2) vs. £ 3 / 2 plane de- 
rived from the total magnetization J z of the ground state. 
Solid lines: for the f 2 -f 3 model. Dashed lines: with no re- 
striction on the occupancy. 

A special case is £ x / 2 = £3/2 = ^5/2 = t (isotropic case). 
Here the ground state is always six-fold degenerate with 
J = 5/2. This is plausible in the limit of large values 
of £: all electrons will occupy orbitals with different j z 
in order to maximize the kinetic energy. With a total of 
five electrons one of the j z orbitals will remain empty, 
implying a six- fold degenerate ground state. 

Even when Uj — > 00, implying the f 2 -f 3 model, but 
with a difference |AZ7j| of any pair of Uj values much 
larger than £, i.e.,\AUj\/t — > 0, one can see that J = 5/2. 
In that case the six-fold degenerate ground states can be 
explicitely written in the form 



lim \y jx ) 



\AUj\/t-+0 



n 



c),(a) + c),{b) 
V2 



|0>, (5) 



where Vp.p is the projector allowing only for atomic f 2 
and f 3 configurations. One implication is that the kinetic 
energy is only — 3£ instead of — 5£ as for free electrons, be- 
cause at most three out of five electrons find free places 
to hop onto. Note that j z 7^ j z t in Eq. (5). While for 
finite values of t/\AUj\ the ground-state wavefunctions 
become vastly more complex, they remain eigenfunctions 
of J = 5/2. Simple results are found again in the limit 
£ — > 0, see Appendix A. 
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As one moves away from the isotropic line with J — 
5/2 the degeneracy of different J z states is lifted. De- 
pending on £3/2 and t 1 / 2 = £5/2 one of them becomes the 
ground state while J is no longer preserved (see Fig. 2). 

The phase with J z = 15/2 appears when t 3 / 2 ^> ^1/2 (= 
£5/2)- It corresponds to a very simple eigenstate in which 
all occupied orbit als have j z > 0. One finds 



= 37 ^ 37 4 /2 (a)ct /2 (a)4 /2 (6)ct /2 (b) |0) 

(6) 



with kinetic energy — t 3 / 2 . The j z =5/2 and j z =l/2 or- 
bitals are occupied on both sites while = 3/2 is occu- 
pied with probability 1/2 on each of the two sites. Note 
that 1^)15/2 is separately an eigenstate of the Coulomb as 
well as the kinetic energy terms. While the Coulomb en- 
ergy is at its absolute minimum, there are no energy con- 
tributions resulting from t 1 / 2 and t 3 / 2 due to the limita- 
tion to f 2 and f 3 configurations. This can be considered 
as an extreme case of partial localization: the j z =5/2 and 
j z =l/2 states are fully localized and the j z =3/2 state is 
maximally delocalized. When the hopping matrix ele- 
ments increase, a discontinuous jump takes place into 
a state with J z = 5/2. The corresponding increase in 
Coulomb energy is seen in Fig. lb. 

The variation of the ground-state energy along the line 
ti/2 — ^5/2 = can be qualitatively understood from the 
following argument: For sufficiently small hopping ma- 
trix elements t 3 / 2 the high spin J z = 15/2 state is low- 
est in energy because it minimizes the Coulomb energy. 
With increasing values of t 3 / 2 it becomes more favorable 
to partially break up the local Hund correlations in or- 
der to gain kinetic energy by partially populating both 
j z = 3/2 and the j z = —3/2 orbit als. In the limit of large 
hopping, their angular momenta components cancel and 
J z is given by the j z components of the remaining elec- 
trons, i.e., by (§ + \) ± § resulting in J z = 11/2 or 1/2. 
For intermediate hopping the total angular momentum of 
states were both ±3/2 and ±1/2 are occupied is given by 
the remaining 5/2 orbitals. For more details on the oc- 
cupation of the individual j z orbitals in the ground state 
see Appendix B. 




FIG. 3. Total momentum as denned by Eq. (7) for sites 
with local f 3 (upper panel) and f 2 configurations (lower 
panel), respectively. Note the scale and the smallness of the 
deviation from the values 24.75 = ((9/2) + 1) • (9/2) and 
20 = (4 + 1) • 4 corresponding to fully polarized f 3 and f 2 
configurations. 

While the total angular momentum J depends 
strongly on the anisotropy, the local angular momenta 
change very little. Fig. 3 shows for different anisotropics 
the local spins defined in terms of projectors onto config- 
urations with given local /-occupation as 

^=2(3) = E < J N 2 ^=2(3) (n)> . (7) 

n—a,b 



The projector P M selects atomic configurations with \i 
electrons. The small range of Jj^ variation in Fig. 3 is 
remarkable. Hopping elements up to 0.5 eV reduce the 
local moments by less than 5% relative to their maximal 
values 24.75 and 20 for f 2 and f 3 configurations, respec- 
tively. One observes a reduction of the local spin expec- 
tation values at large hopping for both the f 2 and the f 3 
site. A decrease of J e ff,/x=2(3) implies for Uj=4 < £7j=o,2, 
an increase of the Coulomb energy. Thus, Figure 3 
demonstrates the competition between Coulomb inter- 
action and hopping. 
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C. Density of states 



D. Field dependence 



In the absence of any hopping, i.e., 



for 



^3/2=^1/2=^5/2 = 0? the ground state is 180-fold de- 
generate (see Appendix A). At finite isotropic hopping 
this degeneracy is only partially lifted. Then J(a) and 
J(b) are no longer separately conserved, but the total J 
still is. 




FIG. 4. Energy levels and their degeneracies for isotropic 
hopping, tx/2 = £5/2 — £3/2 =t = 0.1 eV. Inset: Anisotropic 
hopping A£ = £3/2 — £1/2 lifts the high degeneracies. 



Thus, along the isotropic line, spherical symmetry still 
leads to rather highly degenerate states. As seen in Fig. 4, 
these states are grouped by Hund's rule coupling, i.e., by 
the energy difference of the values Uj. The different clus- 
ters of energy levels are separated by well-defined gaps 
of about leV. 

When spherical symmetry is broken the J z degeneracy 
is lifted. This is seen from the inset of Fig. 4. It shows 
the evolution of the lowest 66 eigenstates as a function 
of the symmetry breaking parameter At = £3/2 — £1/2- 
The value of £ 3 / 2 + £1/2 = 0.2 is thereby kept fixed. Five 
different groups of levels are noticed which for At = 
correspond to the values (ordered by increasing energy) 
J z = 5/2, 3/2, 7/2, 13/2 and a 34 fold degenerate 
level consisting of J z — 15/2 and 17/2 states. When 
At ^ the J z degenerate states split into doublets. For 
At ~ 0.06 one state from the originally 34-fold degen- 
erate states becomes lowest in energy. This marks the 
transition to the "high-spin" state with J z = 15/2. 



In the following we want to discuss the effect of an in- 
ternal or molecular field on the phase diagram. Fig. 5 
shows the changes in the phase diagram caused by a 
ferromagnetic field of three different magnitudes. One 
notices the appearance of new "high-spin" as well as "in- 
termediate spin" phases [6] . Now the phases with low J z 
values require larger hopping matrix elements in order 
to be realized. The regions of "high-spin" states grow 
with increasing field and for small values of £ 3 / 2 a new 
phase with J z — 17/2 appears. This state is similar to 
the J z = 15/2 state and has the form 
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FIG. 5. Phase diagram of the two-site cluster in the pres- 
ence of a magnetic field: from top to bottom h =0.005, 0.01, 
and 0.05 eV, respectively. 



IV. THREE SITES 



In this section, we extend our analysis to a 3-site clus- 
ter. In order to approach the valency f 2 - 7 present in 
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UPt3, we choose a cluster with eight electrons (average 
/ valency 2.67). This choice also avoids the Kramers 
degeneracy which was present in the 2-site cluster and 
which would again be present in the 3-site model with 7 
electrons, i. e., the other possible electron number with 
mixed valency. We put the sites on a ring order to work 
with a model in which all sites are equivalent. As a con- 
sequence, AFM correlations are subject to frustration. 
Furthermore the sign of the matrix elements tj z 's is now 
relevant: the eigenstates of the kinetic energy term in the 
Hamiltonian consist of a doublet with energy tj z and a 
singlet state with energy —2tj z . Thus, for tj z < oc- 
cupancy of j z orbitals two of the three sites is favored 
whereas for tj z > a single occupancy is. Below, we 
consider both negative and positive hopping parameters, 
assuming the same sign for all tj z . All calculations are 
done in the subspace of only atomic f 2 and f 3 configu- 
rations. 
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FIG. 6. Phase diagram for three sites for negative (left 
panel) and positive hopping (right panel). 



2. The case t jz > 0. 



1. The case tj z < 0. 



The left side of Fig. 6 shows the value of the total 
angular momentum J z of the three-site cluster as a func- 
tion of negative t 1 / 2 = £5/2 an d £3/2- As for the two-site 
cluster, a ferromagnetic high-spin phase appears in the 
region £1/2/^3/2 ^ 1- This phase has J z = 12 and can 
be understood as a straightforward generalization of the 
J z = 15/2 phase of the two-site model. It is an eigen- 
state of both the Coulomb and the kinetic terms. Only 
orbitals with j z > are occupied. The j z = +5/2 and 
+ 1/2 orbitals are filled with one electron per site and the 
remaining two electrons occupy j z = 3/2 orbitals. In con- 
trast to this simple ferromagnetic state, all other phases 
in Fig. 6a have a more complicated structure. Along the 
isotropic line of the phase diagram, tj z = t, the ground 
state has J = 6 for all finite values of t v . 

The phases differ also in their magnetic correla- 
tions. Two of them have well pronounced correlations. 
The J z — 1 phase has antiferromagnetic correlations 
E(nm) (J( n ) ' J( m )) < 0)> whereas the J z = 7 and 
12 phases have ferromagnetic correlations. Correlations 
within the J z — 6 phase are by about a factor 5 smaller 
than those of the J z = 12 phase. The other two phases, 
J z = 4,6 have generally small magnetic correlations, 
which can be either ferro- or antiferromagnetic, depend- 
ing on t 1/2 and t 3 / 2 . 



We now turn to the case with positive hopping matrix 
elements. The ground state along the isotropic line has 
J = 4 for all finite values of t v . Away from the isotropic 
line one finds J z = for t\i 2 = £5/2 > ^3/2 an d J z — 3 
for t 1/2 = t 5/2 < t 3/2 , see Fig. 6b. 

As the kinetic energy penalizes states in which two or 
more electrons on different sites occupy the same j z or- 
bital, the ferromagnetic state which was found for tj z < 
is no longer favorable. Most of the phase diagram is 
dominated by AFM correlations. However, for large 
anisotropics £3/2/^1/2 ^ 1 these correlations weaken con- 
tinously with increasing ratio t 3 i 2 jt\i 2 and vanish near 

h/2 = o. 

V. PARTIAL LOCALIZATION 

In the following we want to demonstrate that strong 
on-site correlations result in a considerable enhancement 
of anisotropics in the bare hopping matrix elements. This 
can lead to a a localization of electrons in orbitals with 
relatively weak hybridization. The latter is effectively 
reduced to zero in those cases. 

In order to quantify the degree of localization or, alter- 
natively, of the reduction of hopping of a given j z orbital 
by local correlations, we calculate the ratio of the j z — 
projected kinetic energy Tj z and the bare matrix element 




2 <*gsi( c ±i» c ±^(™) + ^ c -)i*g S > • 

(nra) ,± 

The ground-state wavefunction |^g S ) contains the strong 
on-site correlations. A small ratio of Tj z /tj z indicates 
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partial suppression of hopping for electrons in the ±j z 
orbitals. Two kinds of correlations may contribute to 
that process. The first one is based on the reduction of 
charge fluctuations to atomic f 2 and f 3 configurations. 
This is a result for large values of Uj and can be stud- 
ied by setting all Uj equal to a value much larger that 
the different tj z . The second one is due to differences in 
the Uj values, i. e. , Uj=4 < Uj=2 < Uj=o- The dif- 
ferences in the Uj values are the basis of Hund's rules. 
Hopping counteracts Hund's rule correlations and vice 
versa. What we want to stress is the fact that those cor- 
relations can lead to a complete suppression of hopping 
channels except for the dominant one which is influenced 
only little. 



A. Two-sites cluster 

Results for the ratios Tj z /tj z are shown in Fig. 7. For 
a better understanding one should consider that figure in 
parallel with Fig. 2. The parameter range in Fig. 7 corre- 
sponds to a straight line in Fig. 2 intersecting the abscissa 
at £3/2 = 0.1 eV and the ordinate at t 1 / 2 = £5/2 = 0.1 eV. 
This line crosses three different phases with J z = 15/2, 
5/2 and 3/2, labeled I, II and III respectively. One no- 
tices that in region I only the dominant hybridization of 
the j z = 1/2 and j z = 5/2 orbitals is completely sup- 
pressed. In regions II and III the correlation effects on 
different orbitals are not so extremely different. In a 
forthcoming paper we show how the effects of the cor- 
relations on the ratios Tj z /tj z are changed when all Uj 
values are assumed to be equal. 

To understand the suppression of the Tj z in the limit 
of very small tj z it is sufficient to treat the hopping per- 
turbatively, this is discussed in Appendix A. 
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B. Three sites 



In the case of a three-sites cluster the sign of the matrix 
elements tj z matters because of the respective degeneracy 
of the ground state. We begin with tj z < 0. 

Shown in Fig. 8 are results for different parameter 
choices. Again, for a better understanding, one should 
consider in parallel the phase diagram in the left hand 
side of Fig. 6. Again, three different phases are involved 
and they are labeled as before by I, II and III in Fig. 8. 
Altogether, the figure looks similar to Fig. 7. There 
is again a regime of total suppression of the t 1 / 2 and 
£5/2 matrix elements (region I) and also in region III the 
smaller of the tj z , i. e., t 3 / 2 is strongly reduced to an 
effective T 3 / 2 . The case of equal Coulomb integrals Uj, 
i. e., AUj = is discussed in an extended version of this 
paper. 

Next we turn to the case of positive tj z > 0. Corre- 
sponding results are shown in Fig. 9 which should be con- 
sidered in connection with the right hand side of Fig. 6. 
There it is seen that two different phases with J z — 3 
and are involved with the chosen parameters tj z . In 
distinction to Figs. 7 and 8 there is no regime of total 
suppression of hopping matrix elements. But a strong 
enhancement of anisotropics by the correlations is no- 
ticeable also here. Again, the overall features are very 
similar to those in Figs. 7 and 8 except that the inter- 
mediate regime is missing here. This is of little surprise, 
since we are dealing with local correlations so that the 
effects should be insensitive to the precise form of the 
chemical environment. 



•la- 



'l/2= - 1 




'3/2=° 
t 5/2=-°- 1 



FIG. 7. Values Tj z /tj z for the two-site cluster along the 
line connecting linearly the points written below the figure. 
Phases with J z = 15/2, 5/2 and 3/2, are labeled with I, II 
and III respectively. 



FIG. 8. Values Tj z /tj z of the three-sites model with tj z < 
along the line connecting linearly the tj z values indicated in 
the figure. 
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FIG. 9. Values Tj z /tj z of the three-sites model with tj z > 
along the line connecting the tj z values written below the fig- 
ures. 

It is instructive to take a closer look at the ground 
state. For £1/2/^3/2 < 1> the ground state is twofold de- 
generate with J z = ±3 (see Fig. 6). When ti/ 2 = £5/2 are 
small, £1/2/^3/2 ^ 1? the f 2 configurations are strongly 
dominated by ones in which the orbitals j z = ±1/2 
and j z = ±5/2 are occupied with a ferromagnetic align- 
ment of the electrons. The admixture of other states 
is less then 1.5 % of the norm for £ 3 / 2 = 0.5 and 
ti/2 — ^5/2 = 0.05. Furthermore, the wavefunction is 
strongly dominated, in fact to 87 %, by configurations 
in which the j z = ±1/2, ±5/2 orbitals form an Ising- 
like frustrated AFM background for mobile j z = ±3/2 
orbitals. We describe this AFM background by states 



I -ltt) := (cl 5/2 (a)cl 1/2 (a)) 
(c\ /2 {b)c\ /2 {b)) (c f 1/2 (c)4 /2 (c)) |0> 



(10) 



(| HI) and I HI) are defined analogously). 

We further introduce operators for itinerant elec- 
trons in the ±3/2 orbitals 



C ] - 

U ±3/2 — 



a± 4 3/2 (a) ± 4 3/2 (6) ± 4 3/2 (c) 



a/2±« 



(11) 



with a number a± of order 1 taking into account that 
due to the background not all three sites are equivalent. 
The ground-state wavefunction (with J z = 3) is then of 
the approximate form 



J z 



3) ^ (|ho; ITT) ± |fto; TIT) ± +|fto; III)) /V3 



where 

|fto;itt> 



f _J>p_p c\_ 3/2 c^_ 3/2 



ITT) • 



(12) 



The projector Vp-p projects onto local f 2 and f 3 con- 
figurations. Furthermore |/&o;TIT)> l^oiTTI) are obtained 
from |/iq; |||) via cyclic permutation of sites. Note that 
the kinetic energy contributions Ti/2 an d T 5 / 2 calculated 
with the wavefunction (12) vanishes exactly. 

As the ground state can be written approximately as a 
sum of simple products of itinerant j z = ±3/2 particles 
and a magnetic background, one expects to find a simi- 
lar form for the lowest excited state. This is indeed the 
case: The first excited state has again J z = ±3 but with 
a fourfold rather than twofold degeneracy. The excited 
states can be very well approximated by a difference of 
any two of the three states |/i ; |||), \h ; |||), |/i ; |||), 
e. g. (|ft ;lTT> - |fto;tlt))>/2. For fixed J z = 3 or 
J z = — 3 these states span a subspace of dimension two. 
We may view the first excited state as an excitation of the 
frustrated AFM background which is dominated by the 
j z = ±5/2, ±1/2 orbitals, while the electrons/holes in the 
j z = ±3/2 orbitals simply adjust to the slow dynamical 
changes of the magnetic background. 

The magnetic correlations of the ground state 
are also interesting. The overall magnetic corre- 
lations in \^j z= s) are weakly ferromagnetic, i. e., 
<*J-,=±3|£<m,n> J(rn)'J(n)\*j z=±3 ) = 3. But because 
of (^=3|E< m ,n> J M • J(n)\* Jz= - S ) = -3 one can 
construct states \^j z= s) ± ^^=3) which have no mag- 
netic correlations. Thus, we have a situation with a local- 
ization mechanism akin to a Mott transition, but with- 
out the magnetic correlations which are usually present 
in Mott systems. 



VI. DISCUSSIONS AND SUMMARY. 

In this investigation we have studied the effects of 
strong intra-atomic correlations on anisotropics in the hy- 
bridization of different 5/ orbitals with the environment. 
For that purpose we considered two-site and three-site 
clusters which allow for a complete numerical diagonal- 
ization of the corresponding Hamiltonian. We demon- 
strated that in particular Hund's rule correlations en- 
hanced strongly anisotropics in the hybridizations. For a 
certain range of parameters this may result in a complete 
suppression of the effective hybridization except for the 
largest one, which remains almost unaffected. This pro- 
vides for a microscopic picture of partial localization of 
5/ electrons, a phenomena, observed in a number of ex- 
periments on U compounds. In fact, previous work on the 
de Haas- van Alphen frequencies and anisotropic effective 
masses in UPt3 and UPt2Al3 could explain the experi- 
mental data very well by assuming, that two of the 5/ 
electrons of the U ions remain localized while a third one 
is itinerant [1,2]. The latter is in an orbital with j z = 3/2 
in both cases, and a j-j coupling scheme is used. Previous 
work on the molecule uranocene U(C8H 8 ) 2 gave a hint on 
partial localization [7], but the present work provides a 
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detailed justification for such a feature. By including an 
external field we obtained a rich magnetic phase diagram. 
We plan to investigate its relevance to compound such as 
UGe2. We note that intra-atomic excitations within the 
system of two localized 5/ electrons might provide the 
required bosonic excitations for itinerant electron attrac- 
tion and Cooper-pair formation [8-10]. One may also 
speculate about the application of the present findings 
to multiband Hubbard models. They are used to de- 
scribe correlations in transition metal compounds [11]. 
This problem is left for a future investigation as well as 
the one of a U impurity embedded in a weakly correlated 
chemical environment. Also a molecular-field description 
of the effects found here, e.g., by introducing auxiliary 
fields to model the correlation effects is an interesting 
subject of future research. In an extended version of this 
paper we plan to discuss also the phenomenon of charge 
disproportionation, which may occur in clusters of the 
type investigated here. 



APPENDIX A: PERTURBATION THEORY 



In the limit tj z — > 0, the 180 states of products of lo- 
cal f 2 configurations |4, J' z ) and local f 3 configurations 
\9/2,J z ) are degenerate. The corresponding states are 
direct products of maximal local J states, i.e., three- 
particle states with J(a)(J(b)) = 9/2 and two-particle 
states with J(b)(J(a)) = 4. 

Finite hopping splits the manifold in first order de- 
generate perturbation theory. The transition matrix el- 
ement for hopping with initial product state |/ 3 ; |, J z ) 



\f 2 -A,Jz) and final state \f 2 ;4, J' z - j z ) \ f 



f3.9 j 



-jz) is 



33 
14 



x (f. A ,J>- jz \ Cj Jf.,p>) 

E r 9/2,J z +j z (r 9/2,J' z v 
' 1 J Z °5/2,j z ;4, J z y^5/2,j z ;4,J' z -j z ) 



(Al) 



where the C?'^ are Clebsch-Gordan coefficients. 

In the plane spanned by £ 3 / 2 and £i/ 2 = £5/2? f° r small 
tj z transitions occur near £1/2/^3/2 — 0-35 and at the 
isotropic line, see Fig. 2. These transitions are repro- 
duced with the perturbative treatment. Furthermore, 
the perturbative analysis also describes the degeneracies 
of the ground states correctly. 

Of particular interest is the behavior along the 
isotropic line t 5 / 2 = £3/2 = ^1/2 = t where the lifting 
of degeneracy can be solved in closed form. In this limit 
J' 2 provides an additional good quantum number. We 
therefore consider basis states which are eigenstates of 
J 2 , i.e., |9/2 4JJ Z ) and |4 9/2 JJ Z ). They are obtained 
as linear combinations of the above-mentioned product 
states |/ 3 ; 9/2, J z - J z )\f 2 ;4, J z ) and \f 2 ;4,J z — J z — 



j z )\f 3 ;9/2, J z + j z ). In this basis, the hopping matrix 
elements depend only upon J and can be expressed in 
terms of the 6-j symbols 



t£ 4 JJ Z I X>j» Cj -,(&) I 4 \jj z ) = r(J) 



14 



(A2) 



where 



T (j) = (-iy 10 



459 
^22 

4 j 1 



(A3) 



accounts for the variation with J of the hopping matrix 
element. 

This result allows for a simple interpretation. In the 
j 2 — f 3 model, the total angular momentum J of the 
5-electron states |9/2 4 J J z ) and |4 9/2 J J z ) 
arises from the coupling of three angular momenta. They 
characterize f 2 configurations at the two sites a and 
b with angular momenta 4 and an additional electron 
which can sit at either site. The angular momentum 
j = 5/2 of this electron couples with that of the corre- 
sponding /^-configurations yielding the angular momen- 
tum J(a) = 9/2 or J(b) = 9/2 of the ^-configurations at 
site a or 6, respectively. Technically speaking, the states 
1 9/2 4 J J z ) and |4 9/2 J J z ) arise from applying 
different coupling schemes to the three above-mentional 
angular momenta. Hopping between the two sites con- 
nects the two types of states. As a consequence, the 
hopping matrix elements can be expressed in terms of 
the 6j-symbols which transform between the two differ- 
ent coupling schemes. 

The variation with J of the hopping matrix elements 
is easily evaluated yielding 



r(J) = (1 + 2J) 



183681 - 8424 J - 8344 J 2 + 160 J 3 + 80 J 4 



1064448 



(A4) 



Here J is restricted to half-integer values 1/2 < J < 
17/2. The variation with J of r(J) is displayed in 
Fig. 10. The maximum value for \t(J)\ is obtained for 
J = 5/2 corresponding to a sixfold degenerate ground 
state in the isotropic case. 
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FIG. 10. Variation of r(J') with J in the isotropic case. 



APPENDIX B: ANALYSIS OF GROUND STATES 



At the phase boundaries in Fig. 2 the jumps of the z 
component of the total magnetization of a two-site cluster 
coincide with changes of the orbital-resolved occupation 
number rtj z (m) in the ground-state wavefunction. There 
are shown as gray-scale plots in Fig. 11. They confirm 
the scenario developed in the body of the text and allow 
for some additional observations. The figure corresponds 
to the £3/2 — ti/2 phase diagram shown in Fig. 2 in an 
infinitesimal ferromagnetic field. For weak hopping with 
£3/2 > h/2 = h/2, the j z = +5/2 and the j z = 1/2 or- 
bitals are fully occupied on both sites, whereas j z = +3/2 
is half-occupied, allowing for a maximal kinetic energy 
gain. The corresponding total angular momentum com- 
ponent J z = 2 • § + 2 • \ + § = 15/2. Similarly for large 
£3/2 values the J z = 11/2 phase is formed in which one 
of the j z = 1/2 electrons is promoted to j z = —3/2. 

In all other phases, the average occupations of all 
orbitals are close to one half or zero. Occupations 
near one half could indicate frequent hopping. Such a 
'paramagnetic' interpretation is appropriate only if t v is 
much larger than the energy differences AUj resulting in 
Hund's rule. The three phases to be discussed are best 
classified by looking at the particular j z orbital which 
remains unoccupied (light-colored areas in the figure). 
Obviously, for the case t 3 / 2 < ^1/2 = ^5/2? the loss of 
hopping energy gain is minimal if j z = —3/2 remains 
empty, yielding J z = 3/2. In the case t 3 / 2 >h/2 = h/2, 
either the orbital j z = — 1/2 or j z = —5/2 is unoccu- 
pied, leading to J z = +1/2 and J z = +5/2, respectively. 
It should be kept in mind that the Coulomb interaction 
mixes various j z orbitals. The concept of 'occupying' j z 
states is therefore only approximate. This is reflected in 
small changes of the occupation within each phase. 




1/2 ~ l 5/2 



FIG. 11. Gray-scale plot of orbital occupations as func- 
tion of hopping matrix elements for j z = +|, — | in the 
presence of an infinitesimal ferromagnetic field. Black: fully 
occupied; white: empty. 
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